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ABSTRACT 

'^Gravity  waves  propagating  at  the  surface  of  a  fluid  of  infinite  depth  are 
considered.  The  problem  is  formulated  in  terms  of  a  series  expansion  due  to 
Havelock.  The  series  is  truncated  after  a  finite  number  of  terms  and  the 
unknown  coefficients  are  found  by  collocation.  It  is  shown  that  this  simple 
numerical  procedure  yields  accurate  results  for  waves  of  arbitrary  steepness. 
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SIRNIPICANCE  AND  EXPLANATION 


Over  the  last  15  years  many  efficient  numerical  schemes  have  been 
developed  to  compute  steep  water  waves •  'Riese  schemes  are  often  based  on 
integro-dif ferential  equation  formulations  or  on  collocation  techniques. 

In  this  paper  we  present  a  new  numerical  approach  based  on  an  expansion 
proposed  by  Havelock  in  1919.  Ihis  scheme  is  very  easy  to  implement  and 
yields  highly  accurate  results  for  waves  of  arbitrary  steepness. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MEC,  and  not  with  the  author  of  this  report. 
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STEBP  GRAVITY  WAVES:  HAVELOCK'S  METHOD  REVISITED 
Jean-Marc  Vandan-Broeck* 


1 .  Introduction 

This  paper  deals  with  the  numerical  computation  of  periodic  two-dimensional  gravity 
waves  propagating  at  the  surface  of  a  fluid  of  infinite  depth.  This  problem  was 
considered  before  by  many  investigators.  Most  of  the  existing  numerical  procedures  belong 
to  one  of  two  main  classes. 

In  the  first  class  the  problem  is  formulated  as  an  integro-dif ferential  equation  for 
the  free  surface  profile.  This  equation  is  discretized  and  solved  numerically  by  Hewton's 
method  (see  for  example  Schwartz  and  Vanden-Broeck  ,  Chan  and  Saffman  ,  Vanden-Broeck 
Schwartz'®,  and  Vanden-Broeck'^). 

In  the  second  class  the  solution  is  represented  by  a  Courier  expansion.  The  unknown 
Fourier  coefficients  are  found  analytically  as  series  in  powers  of  a  parameter  equivalent 
to  the  wave  steepness  (Stokes'®,  Schwartz",  Lonquet-Higgins^ ,  Cokelet^)  or  numerically  by 
collocation  (Chan  and  Saffman®,  Rieneckar  and  Fenton'®). 

Numerical  schemes  of  the  second  class  are  usually  inefficient  to  compute  directly 
steep  waves  because  the  Fourier  coefficients  decay  too  slowly  as  the  wave  height 
approaches  its  maximum.  Accurate  solutions  can  however  be  obtained  indirectly  by 
recasting  the  Fourier  expansion  as  Pad(  approximants  (Schwartz",  Lonquet-Higgins^ , 

Cokelet*) .  On  the  other  hand,  very  steep  waves  can  be  calculated  directly  by  using 
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numerical  schemes  of  the  first  class  (Schwartz  and  Vanden-Brosek  ,  Chen  and  Saffman  , 
Vanden-Broeck  and  Schwartz'®). 
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Successful  numerical  procedures  of  the  second  class  have  been  developed  to  compute 
directly  the  highest  vave  (Mlchell^,  Olfe  and  Rottamin^,  Hllllams^^) .  nie  basic  idea  of 
these  numerical  procedures  is  to  represent  the  solution  by  an  expansion  which  takes  into 
account  the  fact  that  the  highest  wave  has  a  sharp  crest  with  a  120*  degree  angle. 

In  this  paper  we  present  a  numerical  scheme  of  the  second  class  Which  enables  us  to 
compute  directly  very  steep  waves-  Our  procedure  follows  closely  the  work  of  Havelock*’ 

A 

and  includes  as  a  particular  case  Olfe  and  Hootman's  scheme. 


2.  Nunerical  Kasults 


He  consiiler  two-dimensional  periodic  waves  of  wavelength  X  and  phase  velocity  C 
propagating  under  the  influence  of  gravity  g  at  the  surface  of  a  fluid  of  infinite 
depth.  We  choose  a  frame  of  reference  in  %ihich  the  waves  are  steady  and  we  introduce 
dimensionless  variables  by  taking  X  as  the  unit  length  and  C  as  the  unit  velocity. 

The  effects  of  compressibility,  viscosity  and  surface  tension  are  neglected. 

We  introduce  cartesian  coordinates  with  the  x-axis  at  the  mean  water  level  and  the  y- 
axis  directed  vertically  upwards.  Gravity  is  acting  in  the  negative  y-dxrection.  Wext  we 
define  the  complex  potential  f  “  *  ♦  and  the  complex  velocity  W  -  u  -  iv.  Here  g 
i.3  the  potential  function,  i|i  the  stream  function,  u  the  x-component  of  the  velocity 
and  V  the  y-component  of  the  velocity.  Without  lose  of  generality  we  choose  i|>  «  0  on 
the  free  surface  and  9  -  0  at  one  crest. 

The  condition  of  constant  pressure  (p  •  0)  on  the  free  surface  can  be  written 

|wl2  +  iLy-i,  ;(,.0  (1) 


where 


2»C 


2 


“  ■  gX  • 

Following  Stokes'^  we  seek  w  as  an  analytic  function  of  f 
plane  i|(  <  0.  This  function  is  periodic  and  tends  to  one  as 


-  ♦  +  i* 
f  ♦ 


(2) 

in  the  lower  half 
Thus  we  have 


W(f  ♦  1)  -  W(f) 
w*1  as  f*-“. 

We  find  it  convenient  to  eliminate  y  from  (1)  by  differentiating  (1) 
Using  the  identity 


(3) 

(4) 

with  respect  to  9 . 


t  i 
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(5) 


we  obtain 


Following  Cokelet'  we  define  the  amplitude  parameter  e‘ 


by  the  relatxon 


(7) 


-  1  -  |w(0)|^|w(1)|2  . 

For  th«  highest  wave  W(0)  -  0  and  e  -  1.  In  general  e  ranges  between  0  and  1- 
The  relations  (3)  and  (4)  show  that  H  can  be  represented  by  the  following 
expansion: 


W(f)  -1+1  b„, 

n-1 


-i2itnf 


(8) 


Because  of  the  symmetry  of  the  wave  about  f  -  0«  the  coefficients  b^  are  real.  They 
have  to  be  found  to  satisfy  (6)  on  -  0.  This  can  be  achieved  approximately  by  using 
the  collocation  procedure  mentioned  in  the  introduction.  Thus  we  truncate  the  series  in 


(B)  after  N  terms  and  we  introduce  the  N  mesh  points 
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(9) 


4N 

Using  (8)  we  obtain  W(ej)  in  terms  the  coefficients  b^j.  Substituting  these  expressions 

into  (6)  we  obtain  N  nonlinear  algebraic  equations  for  the  N  +  1  unknowns 

.2 


U>b^,...,b^.  Another  equation  is  obtained  by  using  (7)  where  is  specified.  This 

sytem  of  N  +  1  equations  is  solved  by  Ne%rt:on'a  method.  Once  the  coefficient  b^^  are 
found)  the  free  surface  profile  can  be  obtained  by  integrating  numerically  (5). 


In  Table  I  we  present  numerical  values  of  u  versus  e  obtained  with  N  =•  60.  For 


comparison  we  also  show  the  accurate  values  of  u  obtained  by  Cokelet'.  Our  values  agree 


with  those  of  Cokelet  to  5  decimal  places  for  e+  <  0.6.  However,  the  accuracy  of  our 


results  decreases  rapidly  as  approaches  1.  This  is  due  to  the  slow  convergence  of 


the  expansion  (8)  as  the  wave  of  maximum  height  is  approached. 

.2 


The  highest  wave,  (i.e.  e^-  1)  is  characterized  by  a  corner  at  the  crest  with  an 

enclosed  angle  of  120*  (Stokes'^,  Amick  et  al . ^ ) .  Therefore 

W(f )  -  f’/^  as  f  +  0  .  (10) 

a  a 

Foilowinq  Michftll  and  Olfe  and  Itootman  we  compute  the  hiqhest  wave  by  replacing  (8)  by 


W(f)  -  (1  -  +  y  . 


(11) 


i 

0 
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i 
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Table  I I 


Values  of  u  for  0.6  <  <  0.99  obtained 

by  using  (8). 


N  -  60 

(2olcelat 

0.6 

1.12229 

1. 12229 

0.8 

1. 17209 

1. 17093 

0.9 

1 . 20068 

1.19014 

0.94 

1.21684 

1.19404 

0.99 

1.253S8 

1.19329 

The  expansion  (11)  satisfies  (10).  Us  truncate  the  expansion  (11)  after  N  -  1  terns  and 
satisfy  (6)  at  the  N  mesh  points  (9).  This  yields  U  equations  for  the  N  unicnouns 
u ,C^ , . . . This  system  was  first  solved  by  Oife  and  Rootman  •  In  particular  they 
found  u  ~  1.93072.  Ue  have  repeated  the  calculation  and  confirmed  this  value. 

The  previous  considerations  suggest  to  combine  the  advantages  of  (8)  and  (11)  by 
representing  the  solution  by  the  expansion 

M(f)  -  (1  -  Be"^^^( ’-'^(1  +  I  d„o"^^’'"*')  .  (t2) 

n»0 

This  expansion  was  first  oroposed  by  HavelocV^.  As  c  ♦  0,  6  «  0  and  (12)  approaches 

(8) .  furthermore  8*1  as  e*1,  so  that  (12)  includes  (11)  as  a  particular  case. 

We  now  truncate  (12)  after  N  -  1  terms  and  satisfy  (6)  at  the  N  mesh  points 

(9) .  Thus  we  obtain  N  equations  for  the  N  ♦  1  unitnowns  8  ,u  .d^ , . . .  ,d^_^.  The  last 

2 

equation  is  given  by  (7)  where  c  is  specified. 

2 

Numerical  values  of  u  vsrsus  e  for  N  •  60,  80  and  120  are  presented  in  Table 
II.  ITie  values  obtained  by  Cokelet^  are  also  shown  in  the  table.  These  results  indicate 
that  the  scheme  converges  as  N  increases.  Furthermore,  the  procedure  yields  values  as 
accurate  as  those  of  Cokelet^  for  values  of  e  close  to  one.  A  comparison  between  the 


--5- 


values  for  M  -  60  in  Tables  I  and  II,  show  clearly  that  the  expansion  (12)  converges 
much  faster  than  the  expansion  ( 8 ) . 


Table  III 

Values  of  u 

for  0.6  <  e2 

<  0.99  obtained  by  using  (12). 

e2 

N  -  60 

M  -  80 

N  -  120 

Cokelet 

0.6 

1.12229 

1.12229 

1.12229 

1.12229 

0.8 

1. 17096 

1.17094 

1. 17093 

1. 17093 

0.9 

1.19007 

1.19025 

1.19019 

1.19014 

0.94 

1.19310 

1.19367 

1.19409 

1.19404 

0.99 

1.19321 

1.19324 

1.19332 

1.19329 

is  worthwhile 

mentioning 

that  Grant^  and 

Schwartr'^  have 

demonstrated  that 

Havelock  expansion  (12)  produces  the  wrong  type  of  singularities  above  the  fluid  (i-e> 
in  <»  >  0).  This  does  not  of  course  invalidate  the  Havelock  expansion.  In  fact,  our 
numerical  results  show  that  this  expansion  is  rapidly  convergent  inside  the  fluid  and  on 


the  free  surface  (i.e.  in  i|i  <  0) 
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